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ABSTRACT 

In this paper the title problem is studied by using Reissner's transverse 
shear theory. The main purpose of the paper is to Investigate the effect of 
stress-free boundaries on the stress intensity factors in plates under bending. 
Among the results found, particularly interesting are those relating to the 
limiting cases of the crack geometries. The numerical results are given for 
a single internal crack, two collinear cracks, and two edge cracks. Also 
studied is the effect of Poisson's ratio on the stress intensity factors. 


1 . Introduction 

In many relatively thin-walled plate and shell structures through cracks 
may develop as a result of cyclic loading. To analyze this fatigue crack 
propagation process the stress intensity factor calculated from the elastic 
analysis of the structure appears to be the most widely used correlation para- 
meter representing the severity of part-flaw geometry and the intensity of 
applied loads. In plates containing through cracks and subjected to membrane 
loading only, usually the solution obtained by ignoring local three-dimensional 
effects and by assuming the validity of conditions of the generalized plane 
stress seems to be quite adequate. Partly because of the practical importance 
of the problem of plates under membrane loading and partly because of the 
relative simplicity of the related elasticity problems, the two-dimensional 
crack problems have been studied very extensively. Even though in many appli- 
cations the bending components of the external loads are also present, as in, 
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v 1 Permanent address: Dept- of Civil Engineering, Istanbul Technical 

University, Istanbul, Turkey. 


- 1 - 


/{fZ'Z/39? 



for example, transversely loaded plates and structures undergoing flow-induced 
vibrations, the solution of the plate bending problem seems to have been 
carried out only for an infinite plate [l]-[5]. These studies have demonstra- 
ted the importance of transverse shear effects on the stress intensity factors 
and have shown that the bending results are sufficiently different from the 
plane stress results. It is, therefore, worthwhile to investigate the influ- 
ence of finite in-plane dimensions, particularly that of stress-free edges 
on the stress intensity factors in plates undergoing bending. 

The problem considered in this paper is a relatively long rectangular 
plate containing coll inear cracks perpendicular to its long sides. Of parti- 
cular interest is the investigation of the edge cracks and crack-free boun- 
dary interaction. As in [l]-[4] the external loads are assumed to by symmetric 
with respect to the plane of the crack and a transverse shear theory [6], [7] 
is used to formulate the problem. 

2. The Formulation of Bending Problem 

Consider a relatively long flat plate of finite width which contains sym- 
metrically located collinear cracks perpendicular to its sides (Figure 1). 

It is assumed that x 2 = 0 is a plane of symmetry with respect to loading and 
geometry and the problem in the absence of cracks has been solved under the 
given applied loads. Thus, through a proper superposition the crack problem 
may be reduced to a stress perturbation problem in which the self-equil ibrating 
crack surface tractions are the only external loads. Also, it is assumed that 
the plate is acted upon by a sufficiently large tensile membrane load so that 
there is no crack surface interference (on the compression side) in the bend- 
ing problem. Thus, the results given in this paper should be considered 
together with the solution given in [8] where the corresponding generalized 
plane stress problem was studied for the same crack geometry as Figure 1. 

By using the Reissner's transverse shear theory, the basic equations for 
elastic plates under bending may be expressed as follows (see, for example, 

[9] for the general case): 

v**w = 0 » (1) 
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< v 2 n - n = 0 , 

(2) 

tcV 2 i|)-i|)-W = 0, 

(3) 

p = + izy. * 20 - li _ liy an 

e x ax 2 ay ' p y ay 2 ax ’ 

(4) 

M xx - iuT [ 5$ + v 3$ + 2 (1 ' v)2 axay 3 ’ 

(5) 

M yy ~ hA^ kft + v a^ “ 2 (1 ’ v) axay ] ’ - 

(6) 

u _ a (1-v) ro 3 2 \p , k it _.w3 2 n a 2 n\n 
M xy “ ~2hX 7r " L [2 a^ay + 2 (1 " v}( aF “ a^ )] ’ - 

(7) 

V = M + < lo +. It 

v x ax 2 ' ' ay ax * 

(8) 

_ aw < n \ an , m 
v y‘ ay'? (1 ' v) w + # * 

(9) 


The dimensionless quantities which appear in (l)-(9) are defined in Appendix 

/ 

A. The dimensions are given in Figure 1. In the usual notation M-j and V./ 
(i,j = 1,2) are the bending^ and the transverse shear resultants ^ and 0 2 ' 
are the components of the rotation vector, u-j, u 2 , and u^ are the components 
of the displacement vector, and a* is a length parameter representing the 
crack size (a* = a for c > 0, d £ b, a* = d for c = 0, d < b. Figure 1). 

As in the corresponding plane stress problem [8], here it is assumed 
that x-j = 0 is a plane symmetry. Thus, in the perturbation problem under con- 
sideration the solution of the differential equations (l)-(3) may be expressed 
as 

co 

w(x,y) = ~ | ( A-j + yA 2 )e” ay cosax da 

0 
co 

+ j | (C-jCOshgx + C 2 x sinhBx)cosey de , 
o 


( 10 ) 



00 CO 

n(x,y) s ~ j B-je r l y sinox da + ^ j B 2 sinh r 2 x singy dg 


01 ) 


aa 

^(x,y) = § | [-A-j+(2Ka - y)A 2 ]e ay cosax da 


uu 

“■ | [-(C-j+2KgC 2 )coshgx - C 2 x sinhgx]cosgy dg , 


where 


r i = [ “ 2 - Tthf > r 2 - i> 2 + 


( 12 ) 


(13) 


and the unknowns A^ , A 2 and B-j are functions of a, and C-j , C 2 and B 2 are 
functions of g. 

By substituting from (10)-(12) into (4)-( 9) the components of rotation, 
the moment resultants and transverse shear resultants are found to be 

oo 

8 x (x,y) ~ ~ j (-[-Aj+(2ica-y)A 2 ]e ay - B^ e ^Jsinax da 


+ — 
i: 


[-(C 1 +2icgC 2 )gsinhBx - C 2 xgcoshgx - C 2 sinhgx 


o 

- B 2 sinh r 2 x]cosgy dg 


(14) 


CO 

g y (x,y) = “ | [-(-A-j+(2Ka-y)A 2 )ae" ay - ^ B-j 


e ^jcosax da 


+ 2 f 

ir 


[[(C-j+2icgC 2 )coshgx + C 2 x sinhgxjg 


~ ~T B 2 cosfl r 2 x ^ slnB y 


(15) 
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M xx (x>y) 


M yy (x,y) 


M xy (x,y) 


CO 

= jjpr ^ {J ( C( 1 —v)ct 2 A-j - (a 2 (l -v)(2<a-y) - 2va)A 2 ]e 
0 

00 

- r i B ] e r ^)cosaxda + j (-^ B r^cosh r 2 x 

o 

- (l-v)B 2 CjCOshBx - [(1-v)x$ 2 sinhBx 

+ 2B(l+(l-v)K8 2 )coshBx]C 2 )cosey dB}» 

00 

= \ {|([-(l-v)a 2 A^ + (a 2 (l-v)(2ica-y) + 2a)A 2 ]e 

o 

+ ^ ^ V ^ ar i B -| e r ^)cosaX da 


(16) 


oo 

+ (- 8r 2 B 2 cosh r 2 x + B 2 (l-v)C-|COshBx 


+ [(^pr " 2vB)coshBx + (1 -v)b 2 x sinhBx]C 2 )cosBy dB), 


(17) 


_a_ 

hA 


|{|([2o 2 A 1 + (2a 2 y-2a-4Ka 3 )A 2 ]e" ay 
0 

00 

- ^7 (r^ 2 +a 2 )B^e ^Jsinaxda + | [^7 (8 2 +r 2 2 )B 2 sinh r 2 x 


+ 2B 2 C-jSinhBx + (2 b 3 x coshBx + (2B+4icB 3 )sinhBx)C 2 JsinBy dB} 


(18) 
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GO 

V x (x,y) = - | [2ica 2 A2e +^-B^e ^Jsinax da 

0 


oo 

+ ~ | [~z slnh r 2 x + ^"2ic6 2 )C2Sinh6x]cos8y dg , 
0 

(19) 

GO 

Vy(x,y) = - j | [2ica 2 A 2 e “ y + ^- B-je ^jcosax da 
0 


OO 

2 f r 2 

+ - J O B 2 cosh r 2 x + BC-j coshgx 

0 


+ (2i<g 2 coshgx + gsinhgx)C 2 ]singy dg , 

(20) 

where 

-2 - 2 _ 10(a*) 2 

Y k(1-v) h z 

(21) 

Because of the assumed symmetry, it is sufficient to consider the problem 
for 0 £ x-| < b, 0 x 2 < » only. Thus, referring to Figure 1, the boundary 
and symmetry conditions of the problem may be expressed as follows: 

M xx (b’,y) = 0, M xy (b\y) = 0, V x (t),y) = 0, 0 < y < oo , 

(22) 

M xy (0,y) = 0, V x (0,y) = 0, B x (0,y) = 0, 0 < y < =o , 

(23) 

M xy (x,0) = 0, V y (x,0) = 0, 0 < x < b' , 

(24) 

My y (x,0) = m(x) , c' < x < d' , 

(25) 

( x , 0 ) = 0, 0 <_ x < c 1 , d* < x < b 1 , 

(26) 

where the normalized length parameters are defined by 


b' = b/a* , c' = c/a* , d' = d/a* . 

(27) 
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From the expressions (14), (18) and (19) it may be seen that the 
(symmetry) conditions (15) are identically satisfied. By using the five 
homogeneous conditions (22) and (24), five of the unknowns A}, B^, , 

(i = 1,2) may be eliminated. The sixth unknown may then be determined from 
the mixed boundary conditions (25) and (26). By substituting from (18) 
and (20) into the homogeneous conditions (24) we find 

A 1 = 1^ B 1 ’ *2 = - l ? 6 ! ' (28) 

If we now define 


B y (x,+0) = f(x) , 0 < x < b' , 


(29) 


from (15), (26), and (28) it can be shown that (Figure 1) 
d' 

B-j(a) = - 2 | f(t)sinat dt , 


(30) 


and 


d' 


Ai (a) = - f(t)sinatdt , Ag(a) = f f(t)sinat dt. 

c' c 1 


(31) 


By using the expressions (16), (18) and (19) the boundary conditions 
(22) may be reduced to 


\ sinh r 2 ^' “ ^2 2g2|C = 7 


B ! sin b ' a da ’ 


(32) 
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B 2 2 — sinh r 2 b' - 2C-j 0 2 si nheb ' - 2C 2 s[&b , coshBb' 


co 

( l+WsinhBb'] -1ft 


By 2 +2a 2 8 


'{a 2 '+B 2 )(r 1 2 +6 2 ) J °1 


] B, sin b'a da , 


( 33 ) 


B 


2 



cosh r 2 b' 


- Ci(l-v)B 2 coshBb' - C 2 [(l-v)B 2 b'sinhBb* 


♦ 2B(l+(l-v)s 2 K)coshBb'] • - ^ f C- xMfthwi 


a(a 2 -B 2 
' WW1 



COSb'a da . 


(34) 


Equations (30)-(34) indicate that all of the unknowns in the problem can 
be expressed in terms of the new unknown function f(t). It is also seen that 
all of the boundary conditions ( 22 ) -( 26) except (25) are satisfied. The 
equation to determine f(t) may, therefore, be obtained by substituting from 
( 30) -( 34 ) and (17) into (25). From the formulation of the problem one may 
observe that the unknown functions A-j , A 2 and B-| refer to the "infinite" 
plate and should give the kernel found in [4]. Indeed, after some simple 
manipulations it may be shown that 


a*(l-v 2 ) 

2:rhX 4 



4k( 1 -v) 
1+v 


•(t-x) 


1 

(t+x) 


j] 


+ 1+v ^-t-x + t+x ^(irlt+xDlMt 

00 

+ lim j C-|COshBx + (*:S 3 coshBx + y- x sinhsx 
y-»+0 o 

- yy BcoshBx)C 2 - r 2 B B 2 cosh r 2 x]cosBy ds 


= m(x) , c' < x < d' , 


(35) 
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where Kg is the modified Bessel function of the second kind. By solving 
(32)-(34) for Cp Cg and Bg and by substituting into (35) we find 



J_ + _L 

1 -x 1+x 


4<(l-y) 

1+v 


^ ( t-x ) 3 + 


It+xF } 




- k(x,-t)}f(t)dt = m(x), c' < x < d' , 


where the Fredholm kernel k(x,t) is given by 




1+e' 2sx -(2b'-t-x)6 


(36) 


-Zrgx 


2k(1-v) 1+e (a 2 e _ ^ b, " t ^ rz - sr e -( b, - t )3u-( b '-x)r 2 

1+v 13 e 3r 2 e ' e 


-2b'S 


. n 28 2b' 8 2 1+e p w 1|r -2Bx» . 4 r rVl -2BXx 

+ [( TT^ ■ TRT j'; e r2"b'e )(1+e ) + i+^(<e 3 (i+e ) 

+ £ x(l-e- 26x ) - s(l+e' 26x )}] 1 [D 1 e- (2b '- t - x ) B 


+ D 2 e- (b '- x)8 e - (b '- t)r 2] - & B 2r 2 (l-te- 2r 2 x ) 1 [Dje-C 5 '' 1 ) 8 


-2b'B 


♦ o 2 e-( b, -t>>-2] e -( b '-x)r 2 ) de , 


(37) 


>1 * f§ rjO-e- 26 ' 6 , - 20*e‘ 2b ' s ) 

-2b' 8 

+ '^2 [l-(b'-t)Bj-(l-v)[| (b'-t)-KS 2 ](l-e~ 2b ' s ) 


(38) 
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( 39 ) 


„ _ 2g2 l+e" 2b ' r 2 „ „-2b’g x \ ft _-2b% 

0 2 " yZ~ ] — -2 b ' r g (1 ' " Kg ( 1 " v )(‘- e ) » 


D = 4b'6 2 e" 2b ' 6 - 8 + 2K$ 3 )( 1 -e' 4b ' 6 ) 




(40) 


If the cracks are internal cracks as shown in Figure 1, then from (26) 
and (29) it follows that 
d' 

* 

f(t)dt = 0 . (41) 

J . 


Thus, the integral equation (36) must be solved under the additional condition 
(41). 

From the following asymptotic behavior of K 2 (z) for small values of 
z 

k 2 (z) = - \ + 0(z 2 log z) , (42) 

it may be shown that the kernel of the integral equation (36) has only a Cauchy 
type singularity. Hence, the solution is of the form 


f(t) = 


F(t ) 

(t-c'rtd'-tr* 


c 1 < t < d' , 


(43) 


and the bounded function F(t) may be obtained by using the numerical method 
described, for example, in [10]. 

If the plate contains a single symmetrically located crack, i.e., for 
c=0, d<b (see Figures 1 and 3), by using the symmetry of the problem and by 
observing that f(t) = -f(-t), (36), (41), and (43) may be epxressed as 


d' 

d(l-V 2 ) f ,r3+V 1 

2»hx** I U l+V t-x 
-d' 


4<(l-y) 

1+v 




+ 


4 1 

1+v t-x 


K 2 (y|t-x|)] 


+ k(x,t)}f(t)dt = m(x), -d 1 < x < d' , 


(44) 
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( 45 ) 


d' 

| f(t)dt = 0 , 

-d 1 

f(t) = , -d < t < d , (46) 

/d z -t z 

where a* = d is used for the normalizing length parameter (Appendix A). 


3. The Stress Intensity Factors 

In the symmetric plate problem under consideration the bending component 
of the Mode I stress intensity factor at the crack tips is defined by (Figure 1) 

k lc ( x 3 ) = Tim [2(^)1* a 22 (x r 0,x 3 ) , (47) 

-*c 

ki d ( x 3 ) = lim ^(x^d)] 15 a 22 (x 1 ,0,x 3 ) . (48) 

x-j-k! 

Let a b be a stress amplitude calculated on the plate surface and used for 
normalizing the stress intensity factors. For example, in a plate subjected 
to uniform bending = M Q away from the crack region 

°b = 6M o /h2 • (49) 

The stress intensity factors are then normalized with respect to a b > / iT*’. If 
the stress intensity factors on the plate surface are defined by 


k(c) = k lc (h/2) , k(d) = k ld (h/2) , (50) 

it is sufficient to calculate k(c) and k(d) in terms of which we have 

k lc^ x 3^ = h72” k ^ * k ld^ x 3^ = WZ k ^ ’ 
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We now note that (36) gives the normalized bending resultant m(x) on 
y = 0 outside as well as inside the crack. Thus, a relatively straightforward 


asymptotic analysis would show that [11] 

k lc^ x 3^ "WT 1 ™ 3x7’ ( 52 ) 

X*j -*C 1 

k ld^ x 3^ = " W1 T /2 ( d-x i) • (53) 

X-j "K1 I 

From (43) and (51 )-(53) it then follows that 

k < c > = 45^T f ' c '> • k < d > “ - F « d '> • < 54 > 

b b 

In plane elasticity problems for cracks it is known that in the close 
neighborhood of a crack tip x-j = d, x 2 = 0 we have 

tf??(x-| ,0) = k +0(1) , (55) 

“ 1 /2(x r d) 



lim 

x-|-k! 


^(d-x-j) C u 2 ( x i * + 0) “ u 2 (x 1 ,-0)] , 


(56) 


where k-j is the Mode I stress intensity factor, p is the shear modulus, k q = 
4-3v for plane strain and < Q = (3-v)/(l+v) for the generalized plane stress. 
In the symmetric bending problem under consideration the crack surface dis- 
placement is given by 

u 2^ x l ,+ 0» x 3) - x 3^2^ x l ,+ ^) " (57) 


From (48), (53) and (55)-(57) it may be observed that the results found from 
the solution of the plane elasticity and the bending problems given by (53) 
and (55) are identical provided k q is selected as (3-v)/(l+v), (i.e., if the 
value for plane stress rather than for plane strain is used for k q ). Also, 
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as shown in [9] and [11] the transverse shear theory used in the present 
analysis gives an angular distribution for the asymptotic stress state around 
the crack tip which is identical to that found for the plane elasticity 
probl em . 

4. The Edge Cracks 

An important special case of the problem described in Figure 1 is the 
edge cracks for which d=b and c>0. In this case as x and t approach the end 
point d'=b' simultaneously the kernel k(x,t) given by (37) becomes unbounded 
and consequently influences the singular nature of the solution. Since the 
integrand in (37) is bounded in any finite interval in 0 <_ e < °>, the unbounded 
terms in k(x,t) will be due to the asymptotic behavior of the integrand. 

Thus, by separating the asymptotic part of the integrand, (37) may be expressed 
as 

CO oo 

k(x,t) = | [K(x,t,3) - K«o(x,t,e)]de + j Kjx.t.sjde . (58) 

o o 

The first integral in (58) is bounded and the second may be evaluated in 
closed form. After a somewhat lengthy analysis similar to that described in 
[8] we obtain 

oo 

| Kco(x,t,g)d$ = k g (x,t) = 2b' -x-t “ (2b' -x-tj^ + (2b'-x-t)* * 
o 

One may note that k s (x,t) given by (59) is identical to that found for the 
plane edge crack problem given in [8] and, together with the Cauchy kernel 
l/(t-x), constitutes a generalized Cauchy kernel. 

Referring to the definition of f(t) given by (29) and the boundary condi- 
tion (26), it is seen that the condition (41) is not valid for the edge crack 
and, as pointed out in [8], is not needed for the solution of the integral 
equation (36). In this case the generalized Cauchy kernel k g (x,t) = l/(t-x) + 

k (x,t) has the property that k ( x , b ' ) = 0, k (b',t) =0, and f(t) is non- 
s y y 

singular at t=b' . Thus, the numerical solution of the problem is obtained by 
letting 
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5 


(60) 


f(t) = 


F(t) 
/t-c ’ 


c ! < t < b’ , 


5. Results and Discussion 

The problem is solved for three crack geometries shown in Figures 1, 

3 and 4 and for a uniform bending moment M 22 = M o (per unit plate width) 
away from the crack region. The results for a symmetrically located inter- 
nal crack of length 2d are given in Tables 1-3 (see the insert in Figure 3). 
Since the problem has three length parameters, namely h, b and d, the results 
depend on two dimensionless length constants. Table 1 shows the normalized 
stress intensity factor as a function of b/d for fixed values of b/h. One 
may note that, as expected for (b/d) -> 1 the stress intensity factor becomes 
unbounded. Also, for fixed plate dimensions b and h in the other limiting 
case of d -* 0 the stress intensity ratio is seen to approach unity, which 
is the result given by the plane elasticity solution for an infinite medium 
with a through crack. These trends are clearly observed in Figure 2 where 
the asymptotes are indicated by b/h = constant lines. The behavior of the 
solution for the limiting case of (d/h) 0 may also be shown analytically. 

Referring to (44), in this case the problem is one of an infinite plate for 
which the Fredholm kernel k(x,t) is zero. Noting that for d -> 0 y 0, by 
using (42) it may easily be shown that in limit the kernel of the integral 
equation (44) would reduce to l/(t-x). Also, by observing that the crack 
opening displacement on the plate surface is u 2 = Byh/2, m(x) = (^/(f^E) = 
M Q /(h 2 E) = o b /(6E), and X 4 = 12(l-v 2 )d 2 /h 2 , if we replace f(x) by 

f (x) = 36 y /3x = (2d/h)3v/ax , v = u 2 (x-j ,+0,h/2) , (61) 

equation (44) becomes 
d' 

2 t T { t>x at dt = " CT b * “ d ' < x < d ' » ( 62 ) 

-d' 

which is the integral equation for an infinite plane under uniform stress 
and for which the stress intensity factor is k(d) = a./tT , 
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Table 1. Stress intensity factor in a plate of finite width containing 
a symmetrically located single crack and subjected to uniform 
bending (Mq) away from the crack region, v =0,3, a b = 6Mo/h 2 , 
(Fig. 1, c=0) , 



b 

= 4h 

m 

= 6h 

b * 

* 8h 

■a 

= lOh 



b/d 

k(d)/ff b iff 

PM 


b/d 

k(d)/o b /cT 

b/d 

k(d)/a b /T 

-> 00 

+ 1.0 

-> oo 

-> 1.0 

OO 

+ 1.0 

-»* oo 

+ 1,0 

-* oo 

+ 1.0 

40 

0.9887 

40 

0.9680 

60 

0.9676 

80 

0.9676 

100 

0.9675 

20 

0.9697 

20 

0.9231 

30 

0.9218 

40 

0.9213 

50 

0.9210 

10 

0.9296 

10 

0.8530 

15 

0.8491 

20 

0.8476 

25 

0.8469 

5 

0.8747 

6.666 

0.8125 

10 

0.8045 

13.333 

0.8019 

12.5 

0.7716 

2.5 

0.8694 

5 

0.7910 

7.5 

0.7780 

10 

0.7737 

10 

0.7526 

2 

0.9094 

4 

0.7812 

6.0 

0.7620 

8 

0.7556 


0.7332 

1.5 

1.0664 

2 

0.8320 

3.0 

0.7502 

4 

0.7267 


0.7166 

1.04 

3.7426 

1.04 

3.4252 

1.5 

0.9282 

’ 2 

| 

0.7817 


0.7347 

' 

->1.0 

oo 

1 .02 

5.1526 

1.04 

3.2040 

1.04 

3.0449 

2 

0.7702 



1.01 

7.4209 

1.01 

7.1933 

1.01 

6.9928 

1 .04 

2.9256 



->1 . 0 

oo 

+1.0 

-> oo 

+1.0 

oo 

1 .01 

l 

6.8141 



i 






o 

• 

T 

->• oo 


In bending problem the kernel of the integral equation is a function 
of the Poisson's ratio v. Therefore, unlike the plane elasticity problems 
the stress intensity factors in bending are dependent on v. Most of the 
results given in this paper have been calculated for v = 0.3. However, to 
show the influence of v on the stress intensity factors, for a central crack 
and for two edge cracks the results are also given for v = 0, 0.2, and 0.5. 
Table 2 shows the results for the internal crack. It is seen that the stress 
intensity factor slightly increases with increasing Poisson's ratio. 


- 15 - 















Table 2. The effect of Poisson's ratio on the stress intensity factor in 
a plate of finite width containing a single crack and subjected 
to uniform bending, (cT b =6M 0 /h 2 , b/h = 10, Fig. 1, c = 0). 


k(d)/a b /5T 


u/ u 

v = 0 

n 

o 

• 

ro 

v = 0.3 

v = 0.5 

■+ 00 

TITO 

-V 1.0 

->• 1.0 

-* 1.0 

100 

0.9583 

0.9650 

0,9675 

0.9717 

50 

0.9002 

0.9151 

0.9210 

0.9307 

25 

0.8119 

0.8368 

0.8469 

0.8638 

12.5 

0.7278 

0.7587 

0.7716 

0.7938 

10 

0.7074 

0.7392 

0.7526 

0.7758 

7.5 

0.6868 

0.7194 

0.7332 

0.7573 

5 

0.6689 

0.7023 

0.7166 

0.7416 

2.5 

0.6850 

0.7197 

0.7347 

0.7612 

2 

0.7192 

0.7547 

0.7707 

0.7975 

1.04 

2.8152 

2.8916 

2.9258 

2.9881 

1.01 

6.6721 

6.7715 

6.8141 

6.8886 

-* 1.0 

00 

00 

-> CO 

-V co 


Table 3. Stress intensity factor vs. width-to-crack length ratio in a plate 
containing a single crack and subjected to uniform bending, (v=0.3, 
d/h = 1, a b = 6Mo/h 2 . Fig. 1, c=0). 


b/d 

-> 1 2 

4 

8 

10 

20 -*■<*> 

k(d)/0 b /cf 

-v oo 0.9094 

0.7812 

0.7556 

0.7526 

0.7488 0.7475 


Table 3 shows the effect of b/d ratio on k(d) for d/h = 1. Again, as b-xl 
k(d) becomes unbounded, and as b + »> k(d) is seen to approach the infinite 
plate result given in [4]. 

To give some idea about the distribution of the stresses in the plate, 
the bending moment M 22 (x-j,0) = M(x) is given in Figure 3. The result is 
obtained from (44) which shows that the moment is -M Q for 0 <_ x-j < d, and has 
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a singularity at x-| = d + 0, The figure indicates that M 22 is a monotonically 
decreasing function of x^. 

Some results for col linear cracks shown in Figure 1 are given in Table 4. 
One set of results shows the stress intensity factors for fixed crack and 
plate dimensions and for varying crack location. The other set of results 
shows the effect of crack length for a fixed crack location (as determined by 
its midpoint). As c + 0 or 2a/(c+d) -*• 1 it is seen that k(c) becomes unbounded 
which is expected. The somewhat unusual result in this case is the steep rise 
of k(d) to the single central crack value as c + 0. As pointed out in [4], 
a smooth continuation of k(d) as c ->■ 0 would correspond to the "pinched" crack 
solution - the steep rise in k(d) being the result of the relaxation of the 
crack surface rotation 8 2 at x-j = 0 from zero to the single crack value. 

Even though it does not seem to be possible to analyze this phenomenon in the 
bending problem, it can be done for the coll inear cracks in plane elasticity. 

For example, from the expression given for k(d) [12] 

-0}] . k- XZW (63) 

°0 

in an infinite plane containing two collinear cracks along x 2 = 0, c < | x-j | < d 
and subjected to uniform tension a Q away from the crack region, if we consider 
k(d) a function of c, it may be shown that for fixed d 

^1+1 and 4r k(d) - - =° for c + 0 , (64) 

v* 

meaning that for c 0 the approach of k(d) to the single crack value is very 
steep. In (63) K and E are the complete elliptic integrals of first and 
second kind, respectively. 

The results for the edge cracks are shown in Tables 5 and 6 and in Figure 4. 
Table 5 and Figure 4 also show the results obtained from the plane elasticity 
problem for the identical crack geometry. The effect of the Poisson's ratio 
on the stress intensity factor in the plate under bending is shown in Table 6. 
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Table 4, Stress intensity factors in a plate containing two symmetrically 

located collinear cracks and subjected to uniform bending, (Fig, 1, 
v=0,3, a=(d-c)/2, b/h»10, a b =6M Q /h 2 ), 


o 

II 

JZ 

X? 

c+d = b 

Ta 

c+d 

k(c)/a b ^a 

k(d)/a b v^T 

b/a 

k(c)/a b ^a 

k(d)/a b /T 

+1/9 

+1.0170 


-Xx> 

+ 1.0 

+ 1.0 

0.14 

0.7758 

0.7835 

10 

0.76846 

0.76847 

0.15 

0.7720 

0.7765 

8 

0.76058 

0.76060 

0.2 

0.7685 

0.7685 

7.5 

0.75969 

0.75971 

0.3 

0.7764 

0.7720 

6 

0.76299 

0.76292 

0.4 

0.7923 

0.7803 

5 

0.77550 

0.77552 

0.5 

0.8167 

0.7916 

4 

0.81218 

0.81212 

0.6 

0.8519 

0.8058 

3 

0.9442 

0.9432 

0.7 

0.9014 

0.8232 

2.5 

1 .2083 

1.2041 

0.8 

0.9751 

0.8447 

2.08 

3.1343 

3.2639 

0.9 

1.1172 

0.8734 

+ 2 

-y co 

. 

-y oo 

0.95 

1.3018 

0.8944 




0.96 

1.3743 

0.9000 




0.97 

1.4791 

0.9066 




0.98 

1.6516 

0.9148 

I 



0.99 

2.0268 

0.9265 




+1.0 

-+■ co 

+1.0134 





Table 5. Stress intensity factors in a plate of finite width containing sym- 
metric edge cracks and subjected to uniform bending or membrane load- 
ing away from the crack region, a^SN^/h 2 , a m =N /h, v=0.3 (see 
insert in Fig. 3). ‘ ! 


c/b 

Bending: k(c)/ab> / a’ 

Tension 


b=10h 

b=6h 

b=2h 

k ( c ) / cr m *'ct' 

+ 0 

-y GO 

oo 

-y co 

-y oo 

0.01 

0.05 

0.1 

8.7889 

3.4726 

2.2893 

2.4567 

2.7957 

2.9467 

0.2 

1.5754 

1 .6656 

1.9359 

2.1769 

0.3 

1 .3082 

1.3690 

1 .5802 

1 .8744 

0.4 

1.1693 

1.2171 

1.3922 

1.7136 

0.5 

1 .0878 

1.1130 

1.2849 

1.6328 

0.6 

1.0396 

1 .0800 

1.2266 

1 .6080 

0.7 

1.0157 

1.0583 

1.2069 

1 .5970 

0.8 

1.0170 

1 .0665 

1.2320 

1 .5915 

0.9 

1.0694 

1.1369 

1.3501 

1 .5883 

0.95 

1.1666 

1 .2645 

1.4798 


0.98 

1.3466 

1.4383 

1.5452 


0.99 

->1.0 

1.4589 

+1.5869 

+1.5869 

+1.5869 

+1.5869 
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Table 6, Effect of Poisson’s ratio on the stress intensity factors in a plate 
containing symmetric edge cracks which is subjected to uniform bending 
away from the crack region, b/h = 10, a. = 6M /h 2 (Fig. 1, d=b, 

2a = b-c), 00 


c/b 

k(c)/a b ^f 

o 

II 

v = 0.2 

v = 0.3 

v = 0.5 

-> 0 

-> OO 

00 

-V OO 

OO 

0.05 

3.3451 

3.4352 

3.4726 

3.5360 

0.1 

2.1681 

2,2533 

2.2893 

2.3514 

0.2 

1.4744 

1.5452 

1,5754 

1.6281 

0.3 

1.2197 

1.2816 

1.3082 

1 ,3548 

0.4 

1.0889 

1.1450 

1.1693 

1.2118 

0.5 

1.0132 

1 .0653 

1.0878 

1.1274 

0.6 

0.9695 

1 .0184 

1.0396 

1 .0769 

0.7 

0.9495 

0.9956 

1.0157 

1.0510 

0.8 

0.9548 

0,9981 

1.0170 

1.0502 

0.9 

1.0127 

1.0522 

1.0694 

1.1000 

0.95 

1.1159 

1.1151 

1.1666 

1.1943 

0.98 

1 .3003 

1.3293 

1.3466 

1.3631 

- 1 

-> 1 .5869 

+ 1 .5869 

+ 1 .5869 

-> 1.5869 


The results given in the tables are self-explanatory. It should, however, 
be emphasized that (a) the stress intensity factor for the plate under bend- 
ing increases with decreasing b/h ratio, (b) bending values are always smaller 
than those of the plane elasticity, (c) as c ->■ b (or as the crack length 2a 
approaches zero) bending as well as the plane elasticity results approach to 
that of a semi-infinite plane containing an edge crack (i.e., k(c) ■> 1 . 5869a / a \ 
a = a b or a = a ra ) (see [8]), and (d) if the results are normalized with 
respect to a'/c (rather than a/a"), it is seen that as c -> 0 in both cases 
k(c)/a'/£" approach 2/tt which is the value obtained from the closed form elas- 
ticity solution of an infinite plane containing two semi-infinite edge cracks 
and subjected to tension equivalent to an average net section stress a m ', where 

CT m' = CT m b/c ’ a b' = a b b/c ^ see Fl ’9 ure 4 )- 
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APPENDIX A 


The normalized quantities for the plate bending problem, 

* * * 

x = x-,/a , y = x 2 /a , z = x 3 /a ; 

* * * 

u = u-,/a , v = u 2 /a , w = u 3 /a ; 


6 x = B 1 ’ 8 y = 3 2 ; 

Mi I M 22 M, 2 

M xx = FE ’ M yy = PE ’ M xy = PE ; 


CT 11 _ a 22 _ a 12 
a xx - E * a yy “ E ’ a xy E ’ 



B = 


5E 

T2TT+vT 



,u _ 12(l-v 2 )a* 2 

x " ' P“ 


In the problem described by Figure 1 a*=a for 0 < c < d <_ b and a*=d for 
0, d < b. 
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Figure 2. The stress intensity factor in a plate of finite width containing 
a symmetrically located single internal through crack which is 
subjected to uniform bending moment = awa y f rom the crack 
region (see insert in Figure 3); v = 0.3, a h = S\^/h 2 . 



v - 0.3 



Figure 3. Distribution of the bending moment M22 (x-j, 0) = M(x-j) in the plane 
of the crack for a plate containing a single symmetric crack and 
subjected to M22( x l ->0) = -Mg on the crack surface -d < x-j < d. 



1.0 


3.0 



Figure 4. Stress intensity factor in a plate containing two symmetric edge 
cracks which is subjected to uniform bending moment or uniform 
tensile stress a m away from the crack region; v = 0.3, b/h = 10, 
ab = 6M o/h 2 > CT b' = CT b b / c > CT m' = a m b /c where and a m ' are the 
average net section stresses. 
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